m 




KEK-CP-197 



Pion form factor from all-to-all propagators of 
overlap quarks 



JLQCD collaboration: T. Kaneko*"^! H. FukayaS S. Hashimoto"^, H. Matsufuru'^, 
J. NoakI T. OnogI ' and N. Yamada"' ^ 



" High Energy Accelerator Research Organization (KEK), Ibaraki 305-0801, Japan 
* School of High Energy Accelerator Science, The Graduate University for Advanced Studies 
( Sokendai), Ibaraki 305-0801, Japan 

Theoretical Physics Laboratory, RIKEN, Saitama 351-0198, Japan 

Yukawa Institute for Theoretical Physics, Kyoto University, Kyoto 606-8502, Japan 

We report on our calculation of the pion electromagnetic form factor with two-flavors of dynam- 
ical overlap quarks. Gauge configurations are generated using the Iwasaki gauge action on a 
16^ X 32 lattice at the lattice spacing of 0.12 fm with sea quark masses down to ms/6, where Wj 
is the physical strange quark mass. We describe our setup to measure the form factor through 
all-to-all quark propagators and present preliminary results. 



The XXV International Symposium on Lattice Field Theory 
July 30-4 August 2007 
Regensburg, Germany 



* Speaker. 

^E-mail: takashi.kanekoOkek.jp 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http://pos.sissa.it/ 



Pion form factor from all-to-all propagators of overlap quarks 



T. Kaneko 



1. Introduction 

Since the pion plays a central role in low-energy dynamics, understanding its properties is of 
great importance. For the electromagnetic form factor Fj^{q^), precise experimental data are avail- 
able near the zero momentum transfer q^ = 0, where the dependence of Fj^{q^) on the quark mass m 
and can be described by chiral perturbation theory (ChPT) provided that m is sufficiently small. 
A detailed comparison of Fji{q^) on the lattice with ChPT and experiments therefore provides a 
good testing ground for lattice calculations in the chiral regime. An understanding on the applica- 
bility of ChPT to lattice data is also helpful towards a reliable calculation of form factors of K, D 
and B mesons. 

In this article, we report on our on-going calculation of Fji{q^) in two-flavor QCD. We employ 
the overlap fermions, which have the exact chiral symmetry and hence allow us to apply ChPT 
straightforwardly to our chiral extrapolation. The salient feature of this study is that Fji{q^) is 
calculated precisely through all-to-all quark propagators |jl|] for a meaningful comparison with 
ChPT and experiments. 

2. Simulation method 

We simulate QCD with two flavors of degenerate up and down quarks using the Iwasaki gauge 
action and the overlap quark action with the standard Wilson Dirac kernel. To reduce the computa- 
tional cost substantially, (near-)zero modes of the kernel are suppressed by introducing two-flavors 
of unphysical Wilson fermions and twisted mass ghosts [^, which do not change the continuum 
limit. Our numerical simulations are carried out on a A'^^^' x Nt = 16^ x 32 lattice at a single value of 
j8 =2.30. The lattice spacing is a = 0.1 184(16) fm, if ro = 0.49 fm is used as input. We take four 
quark masses m = 0.015, 0.025, 0.035 and 0.050, which cover a range of [m^/6,ms/2]. Our current 
statistics are 50 configurations separated by 100 HMC trajectories at each m. So far, we have simu- 
lated only the trivial topological sector, and effects of the fixed global topology by the extra Wilson 
fermions are to be studied [^]. We refer to Ref.[Q] for further details on our production run. 

3. Measurement through all-to-all propagators 

We construct all-to-all propagators of overlap quarks along the strategy proposed in Ref. [jl|]. 
Low-lying modes of the overlap operator D are determined by the implicitly restarted Lanczos 
algorithm and their contribution to the all-to-all propagator is calculated exactly as 

N 

(D-i)iow = (3.1) 

where ( A ^'^^ u^'^^ ) represents k-th eigenmode and the number of eigenmodes Ngp is set to 100 in this 
study. We note that the overlap operator is normal and the left and right eigenvectors coincide with 
each other. 

The contribution of higher modes is estimated stochastically by the noise method with the 
dilution technique [jl|]. One Z2 noise vector is generated for each configuration, and is diluted into 
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Nd = 3 X 4 X Nt /2 vectors with support on a single value for color and spinor indices and at two 
consecutive time-slices. The high mode contribution 

(0-^)high= Xx^^'r^Wt (3.2) 

d=l 

can be obtained by solving the linear equation for each diluted source 

Dx(^) = (l-Piow)T](^) {d = l,...,Nd), (3.3) 

where d is the index for the dilution and Plow is the projector to the eigenspace spanned by the low 
modes. We employ the four dimensional relaxed CG for our overlap solver [||]. 
In summary, all-to-all quark propagators can be expressed as the matrix 

D ' = £ v« (A^vec = N,p+Nd) (3.4) 

k=i 

constructed from the following two set of vectors v and w: 



w« = {«W,...,m(^-),t](i),...,t](^'')}. (3.5) 



Then, two-point functions with the source (sink) operator at time-slice f^'^ and three-point 
functions with the vector current at t" can be expressed as 



k,l=l 

Crr^vMt"-t,t'-t";P,p')= I <;/(^',p')<'5,(^",p-p')^^r(^'-P)' (3.7) 



k,l,m=l 



where the momentum and smearing function for the initial (final) meson are denoted by p^'^ and 
and 

ff^^^\t,p) = £0(r)w(x + r,O(*^^rv(x,O^')e-'P'' (3.8) 

x,r 

is the meson operator with the Dirac spinor structure F constructed from the v and w vectors. The 
smearing function for the local operator is 0/(r) = 5r,o- 

We prepare the v and w vectors on the IBM BlueGene/L at KEK. The computational cost of 
the determination of low modes is ~ 0.6 TFLOPS • hours per configuration. Solving Eq. ( |3.3| ) is the 
most time-consuming part in our measurement, since it requires Nt/2 times more inversions than 
the conventional method. We observe that, however, the low-mode preconditioning of our overlap 
solver leads to about a factor of 8 speedup and its cost is reduced to ~ 1.7 TFLOPS -hours/conf 
for a given valence quark mass m. The calculation of the meson operator i^^'^\t,p) needs much 
less CPU time than the above two steps: it is about 0.2 GFLOPS - hours/conf for a single choice 
of (m,p,F,0). The calculation of correlation functions are even less costly. These calculations are 
carried out on the Hitachi SRI 1000 and workstations at KEK. 
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The key issue in the all-to-all calculation 
is the re-usability of the all-to-all propagators: 
namely, we do not have to repeat the time-con- 
suming Lanczos step and overlap solver to con- 
struct the meson operator 0'^'^\t,\)) for differ- 
ent choices of (p,r,0)- This is a great advantage 
in studies of form factors, which require an ac- 
curate estimate of relevant correlation functions 
with various choices of the momentum configu- 
ration (p,p')- In this study, we test two smear- 
ing functions <pi{r) and ^^(r) =exp[— 0.4|r|], and 
take 33 choices for the meson momentum p with 
IpI < 2. Note that the lattice momentum is in 
units ofln/L in this article. This setup enables 
us to simulate 11 different values of q^, which 
cover a range of — 1.65 [GeV^] < q^. 

It is also advantageous to average the correlation functions over the momentum configura- 
tions, which give the same value of q^, as well as over the source locations {x,t) with temporal 
separations, namely At = t" — t and At' = t' — t" , kept fixed. This averaging reduces the statistical 
fluctuation remarkably as shown in Fig. 



Figure 1: Jackknife data of three-point function 
'^r5r4r5.0.0..(^'/4,M/4;p,O) with |p| = V2 before 
(top panel) and after averaging over source opera- 
tor locations and momentum configurations (bottom 
panel). Data are normalized by the statistical aver- 
age. 



4. Pion form factor and charge radius 

We calculate effective value of the pion form factor from the ratio 

2Mj, R<j,{At,At';\p\,\p'\,q^] 



F^^^{At,At';q^) 
R^{At,At';\p\,\p'\,q^) 



EM)+EM) (A?, A?'; 0,0,0) ' 



7574r5,<t"l> 



(Af,Af';p,p') 



Crsrs-'lxl'/ (^^; P) C75 75,0,0 (A?';p') ' 



(4.1) 
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Figure 2: Effective value of pion form factor Fx_^{At ,At';q^) at m = 0.025. In the left panels, the data are 
plotted as a function of At' with A? fixed, whereas the right panels show At dependence with At + At' fixed. 
Open (filled) symbols show data with (j)~(pi (0$). 
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Figure 3: Pion form factor at m = 0.015 (left panel) and 0.050 (right panel) as a function of q^. The Solid 
line shows the parametrization of the measured pole plus the quadratic correction. The dashed line is the 
expectation from VMD. 

where = ^/ or 0^, and the pion mass M,,: and energy Eji are determined by single-cosh fits to 
Cyj^j 0,0,. We note that the ratio is calculated from correlation functions averaged over the 
momentum configurations and source locations. 

An example of Fji^^{At,At';q^) is plotted in Fig. The pion form factor Fjt{q^) is determined 
from a constant %lioFT^^^^{At,At';q^) in a range of {At, At'), where Fji^^^ shows a reasonable plateau 
and good agreement with data with ^ = <pi. As shown in Fig. ||, we obtain an accurate estimate of 
Fji[q^) except at our smallest q^, where Cy^y^y^^^^ suffers from the most serious damping factor 
^-£.(|p|)Ar^-£.(|p'|)A/' ^jjj^ (IpI, Ip'I) = (2, 1). 

In the same figure, we observe that the q^ dependence of our data is close to the expectation 
from the vector meson dominance (VMD) hypothesis Fji{q^)r^ 1/(1 — g^/Mp) particularly near 
q^ = 0. The q^ dependence is therefore parametrized by the following form of the vector meson 
pole with a polynomial (up to cubic order) 

1 



or an additional pole correction 



l-q^/Ml 



+ ■ 



l-q^Ml 1 



+ ciq^ +C2q^ + CT, q^, 



{c + c'=\] 



(4.3) 



(4.4) 



pole 



While the simplest form Eq. (4.3) with the linear correction (c2,C3 =0) gives a slightly high value 
of ;^^/dof > 2 at heavier quark masses m > 0.035, other fitting forms describe our data reasonably 
well at all m. 

In Fig. ^ we compare the charge radius 



(r2) = 6dF^{q^)/dq^ 



1^2=0 



(4.5) 



obtained from different fitting forms and ranges for the parametrization of the q^ dependence. 
Our result is quite stable against variation of these fitting setup. In the following, we employ 
Eq. (13) with the quadratic correction, since it gives the least value of ;^^/dof with reasonably 
well-determined fitting parameters. We include the leading finite volume correction into the 
result for (r^). 
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Figure ^ shows our chiral extrapolation 
of (r^). In this preliminary report, we test the 
NLO ChPT formula |§1 

{r^)=co + j^^log[Ml]+c,Ml, (4.6) 

where a higher order analytic correction is in- 
cluded to account for the quark mass depen- 
dence of the contribution of the vector reso- 
nance 6/Mp. With two values of /o from our 
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Figure 4: Charge radius obtained at m ~ 0.025 from 
various choices of fitting form (left panel) and lower 
cut for fit range in parametrization of dependence of 
pKiq^) (right panel). 



Studies in p- and £-regimes Eq. ( |4.6[ ) gives 
a reasonable value of ;^^/dof ~ 1 .2 even with- 
out the higher order term. It is however likely 
that this consistency with NLO ChPT is acci- 
dental, since, as seen Fig. ^, the quark mass dependence of our data is mainly caused by that of the 
resonance contribution. 

This chiral extrapolation leads to our preliminary result 

(r^) = 0.388(9),tat(12)sys (4.7) 

where the systematic error is estimated by changing the parametrization form of the dependence 
of Fj[{q-^) and the choice of /o, and by removing the higher order correction in Eq. (p^. This result 
is significantly smaller than the experimental value 0.452(1 1) fm |^. We need further investigations 
on systematic uncertainties: namely the lattice scale has to be fixed from an experimental input and 
we need study finite volume effects including those due to the fixed topology [^. The consistency 



with ChPT may also be tested within the framework including resonance contributions [IC] as in 



an analysis of experimental data in Ref [11]. 



5. Conclusions 

In this article, we report on our calculation of F^{q^) in two-flavor QCD through all-to-all 
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Figure 5: Left panel: chiral extrapolation of charge radius (r^). The experimental value in Ref.[^ and 
(r^) =6/Mp from VMD are alto plotted. Right panel: comparison of (r^) from recent studies. 
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propagators of the overlap fermions. Our preliminary result for (r^) is significantly smaller than 
experiment as in most of previous studies ^jj pj| , 14, 15, |l6[ [I^ shown in Fig. ^ To understand 
the source of this discrepancy, we are completing our measurement of Fj[{q^) with our full statistics 
(10,000 trajectories at each m) for a more stringent comparison with experiment and ChPT. 

We also observe that the all-to-all propagators provide a very precise determination of meson 
correlators. Our studies are already underway for the pion scalar form factor, K ^ n form factors 
and flavor singlet mesons using meson operators Eq. (3.8) stored on the disk. 



Numerical simulations are performed on Hitachi SRI 1000 and IBM System Blue Gene Solu- 
tion at High Energy Accelerator Research Organization (KEK) under a support of its Large Scale 
Simulation Program (No. 07-16). This work is supported in part by the Grant-in- Aid of the Min- 
istry of Education (No. 17740171, 18034011, 18340075, 18740167, 18840045, 19540286 and 
19740160). 
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